arXiv:1508.04840vl [nucl-th] 20 Aug 2015 


EPJ Web of Conferences will be set by the publisher 
DOI: will be set by the publisher 

© Owned by the authors, published by EDP Sciences, 2015 


The Relativistic Three-Body Bound State in Three-Dimensions 

M. R. Hadizadeh 1 - 3 , Ch. Elster 1 , and W. N. Polyzou 2 

1 1nstitute of Nuclear and Particle Physics and Department of Physics and Astronomy Ohio University, 
Athens, OH 45701, USA 

2 Department of Physics and Astronomy, The University of Iowa, Iowa City, IA 52242, USA 


Abstract. Studying of the relativistic three-body bound state in a three-dimensional 
(3D) approach is a necessary first step in a process to eventually perform scattering cal¬ 
culations at GeV energies, where partial-wave expansions are not useful. To this aim we 
recently studied relativistic effects in the binding energy and for the first time, obtained 
the relativistic 3B wave function jT). The relativistic Faddeev integral equations for the 
bound state are formulated in terms of momentum vectors, and relativistic invariance is 
incorporated within the framework of Poincare invariant quantum mechanics. 


1 Introduction and Formalism 

While most three-nucleon calculations utilize partial-wave (PW) representations, at laboratory ener¬ 
gies above about - 300 MeV those become inefficient. For applications in the GeV regime one needs 
to consider vector formulations to implement the dynamics M. Convergence in the GeV regime 
has been demonstrated m in a relativistic one channel model using this approach. In this work we 
review the formulation of the 3B bound state problem in vector variables, which demonstrates that the 
relativistic problem can be treated using a straightforward generalization of methods that have been 
successfully used in the non-relativistic three-body problem. 

In a Faddeev formulation, the relativistic bound state of three particles with mass m interacting 
with pairwise forces, is described by 

|iA> = (M, - M 0 y l T(M t )P\ifr), (1) 

where M, = E t + 3m is the 3 B mass eigenvalue, Mq is the non-interacting 3 B mass operator, T(z) := 
V + V(z - M) 1 V is the boosted two-body (2B) f-matrix embedded in the 3B Hilbert space, and 
P = PjjP + PfcP is the permutation operator for three identical particles. Projecting Eq. Q on 
relativistic momentum states |p, k) leads to the integral equation, 

<p,k|<A> = J dp' elk' dp" dk" k) ^(p, p';M t - <p / , kjP|p", k") <p", k». (2) 

The relativistic basis states contain two Jacobi momentum vectors p and k. The relativistic Jacobi 
momentum k is constructed by boosting the single particle momentum to the 3B rest frame, whereas 
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the momentum p is obtained by boosting the momentum k to the 2B rest frame. In the kernel of 
relativistic Faddeev integral equation one needs the fully off-shell boosted 2B f-matrix T(z) which 
can be calculated from the non-relativistic f-matrix t m by a two-step process, which guarantees that 
the nonrelativistic and relativistic 2B observables and CM wave functions remain unchanged. Those 
steps are 

1. The relativistic right-half-shell 2B f-matrix is related to the non-relativistic right-half-shell by 
a multiplicative function F from 0, 



(3) 


2. The relativistic fully off-shell f-matrix is obtained from the relativistic right-half-shell f-matrix 
by solving a first resolvent type equation, 



(4) 


Here go(z) is the relativistic two-body propagator, with z being the right-half-shell energy and 
Zj the fully-off-shell energy. This procedure was proposed by Keister and Polyzou 0 and 
implemented for the first time by Lin et al. in relativistic Faddeev equations for 3B scattering 


0 . 


By following this strategy one can directly obtain the fully off-shell boosted f-matrix from nonrela¬ 
tivistic one. 

The permutation operator P can be evaluated in relativistic basis states as 


<p',k'|P|p",k"> = 1 (p , ,k , | P ",k")2+ l< P \k'|p",k">3 

= iV(k',k")j<5 3 (p' -k" - ik'C(k',k"))<i 3 (p" +k' + ik"C(k",k')) 


+ 



(5) 


where 


N( k',k") = AT'f-k' - k", k") AT'f-k' - k", k')- 
Here N is the square root of the Jacobian of the basis change. 


( 6 ) 



^l(k/, k/ ) ^m(Pjk) ^m(P jk) (kj)OJ m (^k) 

f^(P jki k/X') + LJ m (kk) OJ m (pji<)OJm(pjk) 


(7) 


The function C is given by 


C(k, k') = 1 + 


cj,„(k') - a> m ( |k + k'|) 


( 8 ) 



where u> m (k ) = Vm 2 + k 2 . In the large mass limit, i.e. m » |p|, |k|, the permutation coefficient C and 
the Jacobian function N approach the nonrelativistic limit and are equal to one. 
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After evaluating the permutation operator and the matrix elements of the fully off-shell two-body 
f-matrix, the relativistic Faddeev integral equation reads, 

<Kp, k) = —- l—— f dk! N{ k, k') r W;p, 7f; M, - co m {k)\ ifa, k'), (9) 

M, - M 0 (p, k) J V / 

where T sym ( p, p';e) = 7(p, p'; e) + T(-p, p';e) is symmetrized boosted two-body f-matrix and the 
shifted momentum arguments are given by 

n = k' + l -C( k, k) k, n = k + ^C(k', k) k'. (10) 

The relativistic and non-relativistic interaction models are defined so that they have the same two- 
body center-of-momentum scattering cross sections and wave functions. In the three-body space they 
are fixed by S -matrix cluster properties and Poincare or Galilean invariance, respectively. Differences 
between the relativistic and nonrelativistic Faddeev equations arise from: 

1. the Jacobian function N, representing from the change of the 3B basis states, 

2. the coefficient C, which appears in the shifted arguments of the momenta, 

3. the relation between the relativistic and non-relativistic f-matrices defined by function F, 

4. the relations between the relativistic and non-relativistic free Green’s functions. 

2 Results and Discussion 

For our numerical analysis we used two models of a spin-independent Malfliet-Tjon type potential, 
MT-V 0 and MT-Vc @. The model MT-Vc contains an additional cutoff function of dipole type. 
The relativistic Faddeev integral equation is solved using a Lanczos technique |7j. The calculated 
relativistic and nonrelativistic 3B binding energies for MT-V potential are E r = -7.4825 (MeV) and 
E m = -7.7382 (MeV), which indicates a reduction of about (E r - E nr )/E nr = 3.3 % in the 3B binding 
energy. This reduction is consistent with the partial wave calculation of Ref. 0, truncated to s-wave, 
with about 2.70 % reduction in the 3B binding energy. The similar behavior can also be seen for MT- 
Vc potential, where the relativistic effects leads to about 3.2% reduction in 3B binding energy. The 
obtained nonrelativistic and relativistic 3B binding energies for MT-Vc potential are E nr = -7.7177 
(MeV) and E, = -7.4732 (MeV), respectively. 

While our calculations indicate that the overall relativistic effect is small, this results from large 
cancellations in the four relativistic corrections mentioned above. In Fig. 1, we show the contribution 
of each correction to the 3B relativistic binding energy of MT-V potential, where E appr is the 3B 
binding energy when one of the ingredients is replaced by the corresponding non-relativistic quantity. 
By setting the Jacobian function N to one in our relativistic calculations, the 3B binding energy has 
a small increase of about 0.8 %, while setting the permutation coefficient C to one leads to a small 
reduction of about 0.6%. Replacing the relativistic right-half-shell f-matrix with the non-relativistic 
one, i.e. setting F = 1 in Eq. <0 leads to about 2.4 % increase in the 3B binding energy. Finally, 
replacing the relativistic free propagator with the nonrelativistic one, the 3B binding energy is reduced 
by about 1.84%. By setting both Jacobian functions, N and C, to one in our relativistic calculations, 
the net contribution of the dynamic ingredients to the 3B binding energy is a small increase of about 
0.15%. Similarly, the net contribution of the kinematic ingredients obtained by setting F = 1 and 
replacing the relativistic free propagator with the nonrelativistic one, leads to an increase in the 3B 
binding energy of about 0.22 %. 
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Figure 1 . Contributions of 
different relativistic corrections to 
the 3B binding energy of MT-V 
potential. (1) N = 1: the Jacobian 
of the permutation operator (2) 

C = 1 relativistic Jacobi momenta 
—> NR Jacobi momenta, (3) F = 1 
the relativistic half-shell Mnatrix 
—> NR half-shell f-matrix (4) Go: 
free resolvent G 0 —> Gg R . 


The extension of the formulation of the relativistic Faddeev integral equations based on realistic 
nucleon-nucleon (NN) interactions is in progress. The input here is an operator representation of spin- 
dependent NN interaction in a helicity representation, which depends on the total spin and the relative 
momentum in the NN subsystems. Since the 3D approach automatically considers all partial waves 
the number of equations is fixed. The successful implementation in 3 N bound state will pave the way 
to 3 N scattering in the few-GeV energy range. 
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